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Recently  it  has  been  demonstrated  that  three-dimensionality  can  play  an  impor¬ 
tant  role  in  dictating  the  stability  properties  of  any  Gdrtler  vortices  which  a  particular 
boundary  layer  may  support.  According  to  a  linearised  theory  vortices  within  a  high 
Gdrtler  number  flow  can  take  one  of  two  possible  forms  within  a  two-dimensional 
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flow  supplemented  by  a  small  crossflow  of  size  0(i?e  )  where  i?e  K  the  Reynolds 

number  of  the  flow  and  G  the  Gdrtler  number.  Bassom  &  Hall  (1991)  showed  that 
these  forms  are  characterised  by  0(1)  wavenumber  inviscid  disturbances  and  larger. 
0(6'^),  wavenumber  modes  which  are  trapped  within  a  thin  layer  adjacent  to  the 
bounding  surface.  Here  we  concentrate  on  the  latter,  essentially  viscous  vortices  and 
describe  their  weakly  nonlinear  stability  properties  in  the  presence  of  crossflow.  It 
is  shown  conclusively  that  the  effect  of  crossflow  is  to  stcibilise  the  nonlinear  distur 
bailees  and  the  calculations  herein  allow  stable,  finite  amplitude  perturbations  to  be 
found.  Predictions  are  made  concerning  the  likelihood  of  observing  some  of  these  vis¬ 
cous  modes  within  a  practical  setting  and  asymptotic  work  permits  discussion  of  the 
stability  properties  of  modes  with  wavenumbers  which  are  small  relative  to  the  imiilied 
0(G'6  )  scaling. 


*  Research  was  supported  by  the  National  Aeronautics  and  Space  .•\dniiiiistration  under  .NASA  contract 
No.  .NASl-18605  while  the  author  was  in  residence  at  the  Institute  (or  Computer  Applications  in  Science 
and  Engineering  (TCASE),  NASA  Langley  Research  Center,  llam|Uon,  V.A  '23GG5 
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§1  Introduction 


Recently  there  has  been  interest  in  the  effect  of  boundary  layer  growth  on  insta¬ 
bility  mechanisms  to  which  the  boundary  layers  are  susceptible.  Here  concentration  is 
focussed  on  the  Gortler  vortex  mechanism  which  has  been  shown  to  occur  in  both  two- 
and  three-  dimensional  boundary  layer  flows  over  concave  walls.  Much  of  the  early 
theoretical  work  concerned  with  Gortler  vortices  addressed  the  problem  of  describing 
the  linear  stability  of  external  two-dimensional  flows  over  such  boundaries.  Early  con¬ 
tributions  were  made  by,  among  others,  Gortler  (1940),  Smith  (1955)  and  Hainmerlin 
(1956).  Later  Hall  (1982a, b)  argued  that  much  of  this  early  work  is  fundamentally 
flawed  for  all  the  analyses  invoked  the  parallel  flow  approximation  (which  essentittlly 
assumes  that  the  basic  flow  in  which  the  vortices  lie  is  independent  of  the  streamwise 
co-ordinate  and  so  neglects  the  effect  of  boundary  layer  growth).  This  approxinnitiun 
enables  the  linear  stability  equations  to  be  expressed  as  ordinary  differential  equations 
but  Hall  illustrated  that  this  assumption  is  unjustifiable  except  in  the  limit  of  a  small 
vortex  wavelength,  and  indeed  this  is  the  explanation  for  the  considerable  inconsisten¬ 
cies  in  the  results  of  the  previous  studies  (see  the  review  article  Hall  (1990)  for  fuller 
details).  Additionally,  in  the  case  of  a  small  vortex  wavelength  the  Gortler  insttibil- 
ity  may  be  described  by  an  asymptotic  structure  which  takes  account  of  boiitniary 
layer  growth  in  a  rational  manner  and  thence  the  parallel  flow  assumption  is  rendered 
unnecessary  in  the  one  situation  in  which  it  has  any  relevance. 

Hall  (1982a)  examined  boundary  layer  flow  over  the  cylinder  y  =  0,  — cc  <  c  <  oc 
where  the  2-axis  is  a  generator  of  the  cylinder  and  y  measures  the  distance  normal  to 
the  surface.  The  x-  co-ordinate  measures  distance  along  the  curved  surface,  which  is 
taken  to  have  variable  curvature 


{l/b)x{x/l),  (l.lu) 

where  i  is  a  typical  radius  of  curvature  of  the  surface  and  I  is  a  characteristic  lengthscale 
in  the  streamwise  direction.  The  Reynolds  number  Re,  the  curvature  parameter  c^  ;md 
the  Gortler  number  G  are  defined  by 

Re^Uol/v,  6  =  l/h,  G  =  2Ri  6,  (11^) 

where  Uq  is  a  typical  flow  velocity  in  the  x  direction  and  is  the  kinematic  viscosity 
of  the  fluid.  Hall  (1982a)  investiga.ted  the  flow  characteristics  when  Re  is  large  and  c> 
is  small  such  that  in  the  limit  S  — >  0,  G  is  held  fi.xed  at  an  0(1)  value.  By  scaling  the 
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spaiuvise  co-ordinate  r  on  the  honndary  layer  thickness,  it  was  demonstrated  that  fur 
a  non  dimensional  vortex  wavennmber  >  1,  linearised  vortex  modes  within  a  two- 
dimensional  basic  flow  are  neutrally  stable  at  a  Gdrtler  number  G  =  where 

go  is  a  known  0(1)  constant  whose  precise  value  is  dependent  upon  the  properties  of 
the  underlying  basic  flow.  For  larger  wavelengths  the  problem  is  fully  nonparallel  and 
the  linear  stability  equations,  which  now  take  the  form  of  a  set  of  partial  differential 
equations,  have  to  be  solved  numerically,  see  Hall  (19S3).  This  paper  showed  two 
significant  features  of  this  nonparallel  flow  problem,  namely  that  the  ideas  of  a  unique 
stiibility  curve  and  of  unique'  growth  rates  at  a  specified  downstream  location  are 
inapplicable  to  the  Gdrtler  problem  because  the  location  where  a  vortex  commences 
to  grow  is  dependent  upon  the  position  and  the  shape  of  the  imposed  disturbance. 

Hall  (1988)  addressed  questions  concerning  the  development  of  nonlinear  nonpm- 
allel  vortices  within  boundary  layers.  This  numerical  investigation  showed  that  a-^ 
the  nonlinear  disturbance  evolves  the  perturbation  energy  becomes  concentrated  in 
the  fundanu'iital  and  mean  flow  correction;  a  conclusion  consistent  with  the  wetikly 
nonlinear  tln'ory  of  Hall  (19S2b)  valid  for  small  wavelength  vortices.  It  is  well  known 
thiit  Gdrtler  vortices  set  up  by  experimental  means  conserve  their  wavelength  as  they 
move  dowjistream.  Since  flie  botmdary  layer  itself  thickens  it  follows  that  the  loctd 
nondiiiK'nsional  vortex  wavenumber  becomes  large  as  the  vortex  develoi)s.  Thus  the 
small  wavelength  limit  in  the  external  Gdrtler  problem  is  appropriate  to  the  ultimate 
development  of  any  fixed  waveh'ngth  vortex  and  hence  sufficiently  far  downstream  in 
many  flows  the  asymptotic  work  of  Hall  (19S2a,b)  becomes  applicable. 

As  with  all  weakly  nonlinear  investigations  the  results  of  Hall  (19S2b)  are  valid 
only  within  a  neighbourhood  of  the  point  where  the  imposed  perturbation  is  neutrally 
st;il)le.  For  vortices  of  wavenumber  ^  1,  their  development  downstream  of  tlie 
point  of  neutral  stability  is  govenu'd  by  the  solution  of  a  pair  of  coupled  nonlinear 
partial  differential  equations  which  adopt  a  simple  asyminotic  structui'e  at  large  values 
of  A,  where  sX  {X  =  0(1))  denotes  the  distance  of  the  vortex  downstream  of  rhe 
neutral  point.  Formally,  for  large  X .  Hall  A:  Lakin  (1988)  showed  that  this  asymptotic 
structure  could  be  used  to  deduce'  the  flow  configuration  for  fully  nonlinear  Gdrtler 
vortices.  :it  which  point  the  mean  flow  correction  generated  by  the  presence  of  the 
vortices  is  as  large  as  the  basic  (undisturbed)  flow  itself.  The  vortex  structure  derived 
l>y  Hall  X  Lakin  essentially  consists  of  a  core  region  in  which  the  vortex  is  concentrated 
and  which  is  bounded  by  two  thin  layers  in  which  the  vortex  activity  is  reduced  to  icero 
('xponentially.  Further  work  l)y  Hall  &  Seddougui  (1989),  subsequently  reconsidered 
I)y  Bassom  A:  Seddougtii  (1990),  has  shown  that  these  thin  layers  are  susceptible  to 
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secondary  instabilities  which  take  the  form  of  travelling  waves  confined  within  these 
layers.  This  theoretical  investigation,  together  with  the  weakly  nonlinear  account  of 
this  form  of  secondary  instability  described  by  Seddougui  ^  Bassom  (1991)  pro\  ides 
good  qualitative  agreement  with  several  experimental  observations,  notablj'  those  of 
Peerhossaini  (Si  Wesfreid  (198Sa,b). 

All  the  work  described  above  has  addressed  problems  which  arise  wlien  Gbi  t  lei 
vortices  occur  in  two-dimensional  boundary  layers  l>ut  in  many  practical  situations  in 
which  Gortler  vortices  are  known  to  arise  the  basic  boundary  layer  is  three  ■  dimensional. 
For  example,  in  the  case  of  a  boundary  layer  flow  over  a  ol)stacle  or  the  flow  over  a  tur¬ 
bine  blade  the  three-dimensionality  of  the  basic  flow  is  potentially  crucial  and  slioukl 
not  be  neglected.  Most  significantly  the  development  of  laminar  flow  airfoils  has  given 
rise  to  designs  which  have  two  areas  of  concave  curvature  on  the  lower  side  of  the 
airfoil;  when  the  wing  is  swept  the  Iroundary  layer  flow  is  fully  three-dimensional  tmd 
the  previously  mentioned  analyses  are  largely  inapplicable. 

The  first  attempt  to  describe  this  three-dimensionality  effect  was  given  by  Hall 
(1985)  who  examined  the  Gortler  mechanism  in  flow  over  an  infinitely  long  swept  cylin¬ 
der.  The  results  obtained  were  quite  general  and  did  not  require  a  precise  descriiniim 
of  the  particular  boundary  layer  under  investigation.  In  was  shown  in  Hall  (1985  )  that 
it  is  the  relative  size  of  the  crossflow  and  chordwise  flow  over  the  cylinder  which  is  criti¬ 
cal  in  determining  the  vortex  structure.  It  was  demonstrated  that  as  this  rtitio,  stiy 
was  increased  from  zero  the  first  significant  change  in  the  vortex  structure  from  that 
in  the  two-dimensional  case  occurs  when  ~  ^ ,  where  R(,  is  the  (large)  Reynolds 

number  defined  in  (1.1).  In  this  situation  the  vortices  become  time  dependent  tuid 
the  high  wavenumber  modes  no  longer  have  vortex  boundaries  aligned  with  the  iiow 
direction.  Indeed,  as  the  crossflow  increases  further,  the  neutral  vortices  have  ;iX( 
perpendicular  to  the  vortex  lines  of  the  basic  flow.  The  neutral  Gortler  number  f("  the 
vortices  was  predicted  by  a  large  wavenumber  asymptotic  analysis,  the  results  of  which 
suggested  that  for  0(1)  values  of  the  ratio  of  the  crossflow  and  chordwise  ve’  ;_ity  lii  ltls 
the  Gortler  mecharism  is  probably  unimportant  compared  with  Tollmi<  ii-Schlichting 
and  crossflow  type  instabilities.  Indeed,  there  is  some  limited  experimental  cvidimce 
which  supports  this  ’onclusion.  Work  by  Baskaran  Bradshaw  (1988)  has  shown,  at 
least  for  turbulent  boundary  layers,  that  increasing  the  crossflow  velocity  component 
tends  to  destroy  the  Gortler  mechanism. 

Further  investigations  concerning  the  role  played  by  the  crossflow  component  for 
the  Gortler  instability  was  motivated  by  results  obtiiined  by  Denier,  Htdl  k.  Seddougui 
(1991)  (hereafter  referred  to  as  DHS).  The  principal  aim  of  DHS  was  to  provide'  a 
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rational  theory  for  the  receptivity  of  Gortler  vortices  and  specifically  to  explain  how  the 
vortices  may  be  triggered  by  wall  roughness  eleinents.  In  the  course  of  this  work,  DHS 
reconsidered  the  stability  of  a  vortex  in  a  high  Gortler  number  flow  by  implementing  a 
linear,  spatial  stability  analysis.  In  particular,  for  G  ^  they  considered  the  structure 
of  vortices  in  two  different  wavenumber  regimes:  they  reworked  the  analysis  of  Hall 
(1982a)  for  high  wavenumbers  and  also  examined  0(1)  wav'enumber  vortices. 

This  latter  disturbance  mode  is  governed  by  inviscid  equations.  By  considering  the 
region  between  these  two  wavenumber  regimes  DHS  identified  a  new  structure  which 
is  relevant  for  vortices  of  wavenumber  and  which  has  the  property  that  the 

vortices  are  trapped  in  a  thin  layer  of  thickness  0^G“^j  at  the  wall.  Furthermore 

DHS  showed  that  within  this  O^G^^  wavenumber  regime  there  exists  a  unique  must 

unstable  Gortler  vortex  according  to  linear  stability  theory. 

Bassom  ^  Hall  (1991)  (hereafter  referred  to  as  BH)  extended  the  work  of  DHS  to 

consider  the  effect  of  introducing  crossflow  into  the  boundary  layer  flow.  Like  DHS,  BH 

conducted  a  spatial  stability  analysis  of  the  vortex  modes  and  they  demonstrated  that 

for  0(1)  wavenumber  modes  at  large  Gortler  numbers  tlie  crossflow  first  has  significant 

—  1  1 

effects  on  the  two-dimensional  results  once  this  parameter  becomes  0{Re  ^G^).  As 
the  crossflow  increases  the  stationary  vortex  structure  takes  on  an  identity  which  is 
essentially  that  of  a  crossflow  instability;  a  mechanism  first  investigated  by  Gregory, 
Stuart  Sc  Walker  (1955).  BH  surmised  that  for  certain  values  of  crossflow  there  are 
ranges  of  0(1)  wavenumber  space  such  that  vortices  with  wavenumbers  within  these 
bands  cannot  persist.  The  work  contained  in  BH  and  concerned  with  these  essentially 
inviscid  modes  has  been  subsequently  extended  on  two  fronts.  First,  Dando  (1992) 
considered  the  effect  of  allowing  for  compressibility  and  has  shown  that  the  properties 
of  the  modes  are  sensitive  to  the  value  of  the  Mach  number.  Second,  Bassom  (1992) 
investigated  the  role  of  time  dependence  in  the  problem  and  demonstrated  that  al¬ 
though  this  effect  is  important  in  dictating  quantitative  properties,  it  is  not  significant 
in  affecting  the  cpialitative  behaviours  of  the  vortex  structures. 

The  second  part  of  BH  was  concerned  with  obtaining  a  description  of  the  influ¬ 
ence  of  crossflow  on  the  O(G^)  wavenumber  (viscous)  vortices  of  DHS.  It  was  shown 

_  i  3 

that  when  the  ratio  of  crossflow  to  chordwise  flow  becomes  0(/?£  ^G’-)  the  result.s  of 
DHS  need  to  be  modified:  -  significantly  it  was  demonstrated  that  the  introduction  of 
crossflow  into  the  problem  has  a  stabilising  effect,  at  least  according  to  linear  stalril- 
ity  theory.  In  particular,  whereas  in  DHS  it  was  conclusively  proved  tluit  stationary 
vortices  arc  necessarily  tmstable  at  0(G5)  wavenumbers,  this  is  no  longer  true  once 
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crossflow  terms  are  introduced.  In  addition,  for  certain  crossflow  values  there  exist 
neutrally  stable  vortex  modes  whilst  at  large  enough  values  of  the  crossflow  no  vortex 
structures  induced  by  centrifugal  effects  can  exist.  BH  also  investigated  the  stability  of 
linearised  time-dependent  viscous  vortices  and  demonstrated  that  crossflow  also  tends 
to  stabilise  these  modes.  Bassom  (1992)  has  since  shown  how  these  unsteady  viscous 
modes  connect  to  the  corresponding  inviscid  ones. 

The  principal  aim  of  the  present  paper  is  to  extend  BH  to  make  a  study  of  the 
nonlinear  stability  characteristics  of  the  ©(Gs)  wavenumber  modes  elucidated  in  the 
latter  article.  Denier  &  Hall  (1991)  have  made  a  numerical  study  of  the  nonlinear 
evolution  of  these  high  wavenumber  structures  within  a  two-dimensional  boundary 
layer.  They  showed  that  the  nonlinear  equations  governing  the  evolution  of  the  Gortler 
vortex  over  an  0(G~^)  streamwise  lengthscale  are  fully  nonparallel  in  nature.  By 
implementing  the  scheme  described  by  Hall  (1988),  Denier  it  Hall  illustrated  that 
given  a  suitable  initial  perturbation  then  eventually  the  energy  of  the  higher  harmonics 
grows  until  a  singularity  is  encountered  at  some  downstream  position.  It  is  concluded 
that  it  is  this  singularity  which  is  ultimately  the  cause  of  vortices  which  are  originally 
close  to  the  wall  moving  into  the  main  part  of  the  boundary  layer. 

Some  work  has  been  initiated  concerning  the  nonlinear  stability  of  the  inviscid 
modes  within  a  three-dimensional  boundary  layer,  Blackaby  k.  Dando  (private  com¬ 
munication).  When  results  become  available  from  this  work  they  should  prove  inter¬ 
esting. 

In  order  to  make  a  nonlinear  study  of  the  viscous  modes  we  perform  an  analysis 
similar  in  spirit  to  that  performed  by  Seddougui  &  Bassom  (1991)  in  their  investigation 
of  nonlinear  wavy  mode  instabilities  imposed  upon  large,  strongly  nonlinear  Gortler 
structures.  We  have  already  alluded  to  the  fact  that  near  the  right-hand  neutral  branch 
for  Gortler  vortices  Hall  (1982b)  demonstrated  that  the  weakly  nonlinear  vortex  is 
governed  by  a  mean-field  interaction  theory  as  opposed  to  the  more  classical  Stuai  t- 
Watson  (1960)  type  approach.  The  latter  is  much  more  frequently  found  in  nonlinear 
situations  in  a  variety  of  fluid  mechanics  problems.  We  have  noted  that  BH  showed 
that  the  effect  of  increasing  the  crossflow  component  of  the  boundary  layer  flow  is, 
in  general,  to  stabilise  the  vortex  modes.  Indeed  there  are  certain  crossflow/  vortex 
wavenumber  pairings  which  give  rise  to  a  neutrally  stable  structure;  this  structure  is 
confined  to  a  thin  layer  adjacent  to  the  curved  surface  over  which  the  flow  occurs. 
By  gradually  increasing  the  amplitude  of  the  vortices  from  the  infinitesimal  values 
upon  which  linearised  theory  is  based  we  can  identify  the  crucial  vortex  size  at  which 
nonlinearity  becomes  important.  Like  Seddougui  &  Bassom  (1991)  we  find  that  at 


this  point  all  the  harmonics  of  the  disturbance  attain  asymptotically  identical  sizes 
and  then  the  vortex  is  determined  by  the  solution  of  an  infinity  of  coupled  nonlinear 
ordinary  differential  equations. 

There  have  recently  been  a  number  of  problems  studied  which  share  this  charac¬ 
teristic  that  all  the  harmonics  of  the  spamvise  dependence  of  the  vortex  play  important 
roles.  The  resulting  infinity  of  coupled  equations  poses  a  formidable  numerical  chal¬ 
lenge  in  order  to  execute  a  satisfactory  solution  and  to  our  knowledge  there  have  been 
few  attempts  to  tackle  such  a  problem.  One  example  of  such  an  attempt  was  made 
by  Bassom  Blennerhassett  (1992)  who  were  concerned  with  obtaining  the  structure 
of  strongly  nonlinear  vortices  within  curved  channel  flows  and  they  have  shown  that 
at  sufficiently  high  Taylor  numbers  the  flow  structure  develops  a  wall  layer  which  is 
of  the  same  size  as  the  vortex  wavelength.  In  this  case  an  infinite-dimensional  set 
of  differential  equations  is  obtained,  some  properties  of  which  were  investigated  com¬ 
putationally.  Earlier  work  by  Denier  (1992)  concerned  with  the  development  of  fully 
nonlinear  Taylor  vortices  has  given  rise  to  a  similar  set  of  equations  but  the  discussion 
was  restricted  to  explaining  why  these  equations  are  guaranteed  to  have  a  solution 
with  the  desired  boundary  conditions  rather  than  seeking  their  solutions  explicitly. 

Given  that  the  nonlinear  development  of  the  high-wavenumber  ©(G^)  Gortler 
vortex  is  governed  by  the  infinity  of  equations  here  we  address  the  weakly  nonlinear 
limit.  A  Stuart-Watson  type  theory  is  applicable  which  leads  to  an  amplitude  equation 
for  the  evolution  of  the  vortex.  Some  comments  are  made  concerning  the  solutions  of 
this  equation  for  a  variety  of  vort('x  wavenumbers  and  frequencies.  We  note  here  the 
fact  that  for  these  high-wavenuml)er  viscous  modes  the  governing  equations  may  be 
easily  scaled  so  that  the  results  cU’iived  are  valid  for  a  wide  class  of  three-dimensional 
boundary  layers.  This  is  in  contrast  to  the  situation  for  the  0(1)  wavenumber  invis- 
cid  modes  for  which  the  quantit.ifive  behaviour  of  the  vortex  structures  in  entirely 
dependent  on  the  particular  boundary  layer  under  consideration. 

The  reniaind(’r  of  the  paper  is  divided  as  follows.  In  the  coming  section  the 
problem  is  formulated  and  the  ftdly  nonlinear  coupled  set  of  governing  differential 
equations  given.  The  weakly  nonlinear  limit  of  these  equations  is  obtained  in  section  3 
and  the  numerical  methods  iinpleniented  in  their  solution  described  in  section  4.  The 
nonlinear  propertit's  of  the  modes  are  given  in  section  5  and  these  are  then  analysed 
within  a  low  wavenumber  limit.  W<'  close  with  some  discussion. 
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§2  Formulation 

Our  aim  is  to  derive  the  equations  which  determine  nonlinear  high  wavenumber 
viscous  vortex  modes  in  a  slightly  three-dimensional  boundary  layer.  Following  Hall 
(1985)  we  consider  the  boundary  layer  flowing  over  the  cylinder  y  =  0,  —  oc  <  x  < 
oo,  where  the  2-axis  is  a  generator  of  the  cylinder,  y  is  the  distance  normal  to  tlie 
surface  and  the  x  coordinate  measures  distance  along  the  curved  surface.  The  Reynolds 
number  Rg,  the  0(1)  Gortler  number  G  and  the  curvature  parameter  6  are  as  given 
in  (1.1)  where  Uq  is  a  characteristic  flow  velocity  in  the  2--directioii.  The  basic  flow 
within  the  three-dimensional  boundary  layer  is  assumed  to  be  of  the  form 

u  =  Uo{uiX,Y),R:h,{X,Y),R7h*w{X,Y)^  (l  +  o(i?r"))  (2.1) 

I 

where  X  —  x/l  and  Y'  =  yRg  /I,  and  the  crossflow  parameter  A*  is  of  order  one.  This 
boundary  layer  profile  is  in  general  found  by  numerically  integrating  the  downstream 
momentum  and  continuity  equations  to  determine  u  and  e,  and  then  finding  w  from 
the  spanwise  momentum  equation. 

I 

It  is  convenient  to  deflne  the  scaled  spanwise  coordinate  Z  —  Rt  ~/l  and  let  t  be 
the  temporal  variable  scaled  on  I/Uq.  The  basic  velocity  profile  is  perturbed  by 

(u  {t,  X,  Y,Z),  R:  '  H  (f ,  A',  r,  Z) ,  R:  '  H'  (f,  A,  T,  Z))  ,  (2.2u ) 

the  pressure  field  is  given  by 

p  =  p  +  i?7‘P(f,A,r,Z)  (2.26) 


and  on  substituting  (2.2)  into  the  Navier-Stokes  equations  then 


Ux  +  Vy  +  IVz  =  0, 

—  Ut  +  Pyy  +  Gzz  —  uyV  —  TiUx—TixU  —  vUy  —  X*TcU: 

=  UUx  +  VUy  +  WUz, 

—  Vt  +  Vyy  +  \'zz  —  GxuU  —  Py  —  uVx  —  Ty\U  —  uVy  —  cyV'  —  A*ielA 

=  UVx  +  FVV  +  H'Fz  +  ^\UV 

—  \Vt  +  ITyy  +  Wzz  ~  Pz  ~  all  ,Y  —  \*WxU  —  FllV  ~  X*\'wy  —  A*tFll7 

-  UWx  +  VTTy  +  llTlz. 


(2,3a) 

(2.36) 

(2.3e) 

(2.3J) 


Here  terms  of  relative  order  Re  ^  have  been  neglected. 

The  linearised  solutions  studied  by  BH  are  obtained  by  considering  (2.3)  when  all 
the  right  hand  sides  of  (2.36  —  d)  are  put  equal  to  zero.  Additionally,  the  nonlinear 
equations  examined  by  Denier  &  Hall  (1991)  to  determine  the  development  of  large 
vortices  in  two-dimensional  boundary  layers  can  be  retrieved  by  setting  A*  =  0  in  the 
above. 

We  now  follow  the  scalings  first  derived  by  DHS  in  which  it  was  shown  that  the 
O(G^)  wavenumber  vortices  are  confined  to  a  layer  of  thickness  0(G'“s)  adjacent  to 
the  cylinder.  These  modes  have  a  spatial  growth  rate  0{G^),  and  it  was  proved  in 
BH  that  the  three-dimensionality  of  the  basic  flow  is  first  significant  for  these  viscous 
vortices  once  the  scaled  crossflow  A*  is  0(G^ ).  Therefore  it  is  convenient  to  define  the 
0(1)  crossflow  parameter  A  by 

A*=gU.  (2.4a) 

If  the  spanwise  wavenumber  of  the  fundamental  vortex  is  taken  to  be  a  (which,  recall, 
is  0(G®))  then  the  scaled  wavenumber  ko  is  defined  according  to 

a  =  koG'^,  (2.46) 


and  the  disturbance  is  confined  to  the  region  where  the  co-ordinate  ip  given  by 


iP  =  koG^Y, 

is  of  order  one.  In  this  layer  the  nonlinear  disturbance  assumes  the  form 

(n  =  oo  n  =  oo  \ 

^  Uo„El*  +  G-’^  ^  LhnE'P + 

n^—oo  »  =  — oo  / 

(n  =  oc  n=oo 

n=  — oo  n  =  — oo 

(n  =  oo  /J=oo 

Y,  wo^E^  +  G-'^  Y  + 

n=  — oo 


(2.4c) 


w=— OO 

n  =  oo 


p-agU  y  ^on^r  +  G-"  Y  - 


(2.:a) 

(2.56) 

(2.5c) 

(2.5d) 
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where 


£;i=exp  ikoG-^Z  +  G"^  y*^  (^o  +  G"Mi(X)  +  ...)c/X 

-iCt  J  +  +  • 


(2.5e) 


This  definition  of  Ei  refiects  the  fact  that  the  vortex  has  spanwise  wavenumber  0(G^ ) 

2 

and  that  the  results  of  BH  demonstrate  that  the  disturbance  has  an  ©(G®)  growth 
rate  (since  it  transpires  that  jSq  is  purely  imaginary).  Furthermore,  the  parameter  A 
is  at  our  disposal.  Initially  the  choice  A  =  1  is  .made  which  permits  the  fully  nonlinear 
vortex  equations  to  be  deduced  but  subsequentlj'  A  is  allowed  to  tend  to  zero  which 
corresponds  to  taking  the  appropriate  weakly  nonlinear  limit. 

For  F  1  the  basic  flow  quantities  u,  w  are  assumed  to  take  the  forms 


S  =  iin{X)Y  +  + 

w  =  ^,2,(X)Y  +  + 


{■2.6a) 


(2.66) 


and  V  =0(F^).  The  vortex  equations  are  now  derived  by  substituting  (2.4-6)  in 
equations  (2.3)  and  comparing  coefficients  of  like  terms.  The  continuity  equation 
(2.3a)  yields 

0oUon  +  ikoWon  =0,  if  n  7^  0  (2.7a) 


Voo  =  0, 

^  {^lUon  +  0oUin  +  ikoWln)  +  ^’O  ,  ,  =  0. 

dip 


(2.7a) 

(2.76) 

(2.7c) 


The  streamwise  and  spanwise  momentum  equations  {2.3b,  d)  then  give 

^qUoo  +  tT’oWoo  =  —  (--S) 

but  the  constraint  that  the  disturbance  be  confined  to  the  thin  wall  layer  requires  that 
remain  bounded  a,s  ip  oo  which  in  turn  forces  the  relationship 


+  iXfJi2i  —  b. 


(2.9a) 


This  confirms  the  assertion  that  /Sq  is  purely  imaginary  and  then  (2.8)  becomes 


0oUoo  +  ifcoVFoo  —  0. 


(2.96) 
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Next  order  terms  in  the  streamwise  momentum  equation  leads  to  the  equation 


d2 


2  ,  i^oN  iA  (^11/122  - /ii2^2l)  /2  1  rr 


A’o 


oo 


dU^ 

^  dtp  ^ 

n  =  — oo.n^iV  ^ 


O(N-ri) 


n  JJ  dVQ(isi-n) 
e^Ori 


AT  -  /2 


dtp 


+ 


NUo! 


N 


{PoUio  +  /3i  Uqq  +  iAoVKia)  • 


(2.10) 

Our  aim  is  to  derive  the  governing  equations  for  the  leading  order  vortex  compo¬ 
nents  {t'^oni  Aon}-  Clearly  one  relation  between  these  unknowns  is  established  by  (2.10) 
once  an  expression  for  {/SqUio  f^iUoo  -f  iA'oHho)  is  found.  Third  order  components  of 
momentum  equations  (2.36,  d)  may  be  combined  to  show  that  if 


=  /3oU\o  +  iA'olV'^io, 


(2.11a) 


then 


^n  =  —  oo,n^0 

Employing  relation  (2.7c)  and  using  (2.10)  for  iV  =  0  gives 


(2.116) 


{0qU\q  +  idoUoo  -f-iAolEio)  —  —  n  (^  /  dip  ’ 

n=-oo,n5i^0  ”  \o  / 


which  in  turn  reduces  (2.10)  to 


(2.12) 

The  second  equation  for  the  h'ading  order  vortex  components  is  obtained  by  fol¬ 
lowing  the  scheme  outlined  in  BH.  Order  terms  in  the  momentum  equation 
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(2.3c)  enables  a  relation  for  the  leading  order  pressure  gradient  to  be  derived  in  the 
form 


dll)  “  ydil>‘^  2klfin 

-  {f^oUyo  +  PiUoo  +  ikoWio) 

Kq  Kq 


(2.13) 


Also  O(G^)  and  0(G»)  terms  in  the  streamwise  momentum  equation  (2.36)  and  the 
spanwise  momentum  equation  (2.3d)  respectively  show  that 


d4>'^ 


-N^  + 


iQoN  0ii.iuNip  iA(pnp22  -  ^2i/^i2)A^0' 


"^0 


'"0 


2AoV 


11 


^7ia 


finViN 


il’2 


=  RHSu 


(2.14a) 


and 


2  iOfoN 


+ 


/?i/ziiiVt/)  iA(/xn/i22  - /^2i/^i2)-A"i/’^ 


2.3 


2H 


6^11 


=  RHS,, 

Kn  Kq 


(2.146) 


where  the  complicated  terms  RHSi  and  RHS2  involve  {l/on,  Vo,i }.  By  adding  j'o 
multiples  of  (2.14a),  iko  multiples  of  (2.146)  and  multiples  of  (2.12),  differentiating 
the  resulting  equation  with  respect  to  ip  and  substituting  for  dPo,\/dtl’  according  to 
(2.13)  leads  to  the  required  second  equation  relating  the  leading  order  vortex  quantities. 

As  in  BH  it  is  convenient  to  invoke  the  scalings 


ko  =  (xo/2n)'  fc,  d  =  l-o.v  =  ,\o 


koA2  —  \6/^ri^2V»  —  Xudu^i 

A(/2u/^22  -  /221/212) 


2/^fi.\o 


(2.15a  —  c) 
(2.15d-  e) 
(2.15/ -y) 
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which  reduce  the  pair  of  equa^^ions  for  the  vortex  to  the  forms 


f  (P  ^^2  \  jr  ^A' 

Vd^/>2  '  A;3  p  )  ^Ar  ^2 

_  1  ^  fdUn,,  n  dVN-n\  , 

n=-oo,n^N  ^  ^  ^  ' 


where 


1  ^  /iV2(2n-A0dF„  N^dVr,dH^N-n  N  d^Vn\ 

k  ^  [  n  dik  n  dip  dip^  n  dip^  )  ' 

n  =  -oo,n?‘0,iV  A  'r  y  y  / 

(2.16c) 

Notice  that  the  scaling  (2.156)  ensures  that  neutrally  stable  modes  correspond  to  the 
case  when  the  scaled  parameter  j3  is  purely  real  valued. 

These  equations  dictate  the  strongly  nonlinear  evolution  of  vortex  modes  and 
are  the  three-dimensional  extensions  those  equations  solved  by  Denier  &z  Hall  (1991). 
A  numerical  solution  of  the  full  three-dimensional  equations  (2.16)  is  planned  for  a 
future  paper  but  for  present  attention  is  restricted  to  the  consideration  of  the  weakly 
nonlinear  properties  of  these  equations.  Such  a  treatment  enables  predictions  to  be 
made  concerning  the  properties  of  Gdrtler  vortices  in  three-dimensional  bouiidaiy 
layers  of  practical  importance.  It  is  to  be  noted  that  the  scalings  (2.15)  ensure  that 
the  ensuing  deductions  may  be  related  to  a  wide  variety  of  boundary  layers  and  the 
precise  properties  of  vortices  within  any  particular  flow  may  be  easily  deduced  by 
appealing  to  the  recipes  (2.15)  with  the  particular  boundary  layer  functions  xo,  /in, 
k-12,  fJ-21  find  /i22  inserted. 

Consequently,  in  the  following  section  we  concentrate  on  the  case  when  |I7Af|  and 
|F;v|  are  small  or,  equivalently,  in  the  context  of  the  original  expansions  (2.5),  when 
the  parameter  A  -C  1. 
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§3  The  weakly  nonlinear  theory 

Our  analysis  of  the  full  equations  (2.16)  is  concerned  with  the  problem  in  which 
the  vortex  mode  is  of  amplitude  0(/i),  with  h  <C  1,  relative  to  the  scalings  implied  in 
(2.5).  In  this  case,  in  keeping  with  the  usual  weakly  nonlinear  approach  we  anticipate 
from  (2.16)  that  the  perturbation  quantities  are  changed  by  0{h^)  from  their  linear 
values  and  so 


U\  —  hUio  ‘ '  t 

{3.1a, b) 

U2  =  h'^U20  +  ---, 

V2  —  h^V2o  +  •  •  • , 

(3.1c,  d) 

Uq  =  h^Uoo  +  •  •  • , 

$  =  h‘^^0  +  •  ■  • , 

(3.1c,/) 

and  the  scaled  parameters  involved  take  the  forms 

^  —  1^0  +  1^2  +  •  •  • , 

Q  =  Qo  +  +  ■■■■ 

i3.lcj,h) 

The  remaining  terms  in  (2.16)  are  sufficiently  small  so  as  to  become  negligible  in  the 
ensuing  analysis. 

The  expressions  (3.1)  when  inserted  in  (2.16)  suggest  that  it  is  convenient  to  define 
the  operator  hy 

r  _  ilSoNip  iXNip'^ 

=  - P - —■  (3..) 

At  leading  order  in  (2.16)  we  obtain  the  system 

ii  (Cio)  -  ^  =  0.  (3  3«) 

which  needs  to  be  solved  subject  to  the  conditions 

^10,  1^10,  — ^0  as  0 — >0,oo,  (3.3c) 

dip 

in  order  to  ensure  zero  disturbance  quantities  on  the  surface  T  =  0  and  that  the 
vortex  decays  as  Y  — >  oo.  The  system  (3.3)  is,  of  course,  just  the  linear  viscous 
equations  discussed  in  BH.  The  salient  properties  of  the  relevant  solutions  of  (3.3)  will 
be  discussed  in  a  forthcoming  section. 
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where  here,  and  in  the  remainder  of  this  article,  an  asterisk  used  as  a  superscript  on 
a  quantity  denoted  the  complex  conjugate  of  that  quantity.  Also,  at  the  same  order, 
the  harmonic  terms  1/20  and  1-20  are  determined  by  the  pair 


1  I  dUiof.  fj  dViQ 
no  —  i^lO' 


dll’ 


i,  (jL  _  4j 

\d4)'‘  J  ™  A-  \  dip  dp'‘ 


(3.5a) 


_ fr'2 

j^.2  It)’ 


(3.56) 


subject  to  the  boundary  conditions  that 


U20,  V2Q, 


dV, 


20 


dll' 


0  as  V’  — >  0, 00. 


(3.5c) 


The  numerical  procedures  implemented  in  order  to  solve  (3.3)  and  (3.5)  are  de¬ 
scribed  in  the  coming  section  and  so  here  we  concentrate  on  deriving  the  required 
amplitude  equation.  In  order  to  obtain  the  desired  equation  then  it  is  necessary  to 
consider  0{h^)  terms  in  (2.16)  which  lead  to  the  pair 


C.  (Cn)-^=( 


i^2 


Ar2 


A-3 

^1*0  dF2o 
2  dll' 


JT  I  1  I  I  ^^^00  dU20^,, 


+  2U: 


d\\ 


20 


10 


dil> 


-f 


A-2 


(3.6a) 


j  ,  d'  \  .  2,X.  pUn 


1 


/Q2  ,  il324> 


k'^ 


+ 


1-3 


d^ 

dp 


-1  Vio 


1'*'“  (^  ' ')  ~ 


1  rfF.'o  ,  Zf,.dV2„  dV,%SV2„  ,  1  d’VVo 


+  ;v', 


10 


dij’ 


,  I'Vo  dn>3„ 

f/0  0?^)^  '  2  (10  d02  dp 


+  « 


dp  j 

(3.66) 
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subject  to  the  boundary  conditions 

^11)  ^iii  0  as  'ip — >0,oo.  (3.6t) 

dip 

The  homogeneous  form  of  (3.6)  is  merely  (3.3)  and  so,  as  is  the  normal  method 
within  a  weakly  nonlinear  analysis,  equations  (3.6)  only  have  a  suitable  solution  if  a 
certain  compatibility  condition  is  satisfied.  This  condition  leads  to  the  specification 
of  the  correction  terms  ^2  and  within  the  streamwise  dependence  and  frequencj' 
expansions  detailed  by  (3.1^,  h).  To  derive  the  compatibility  requirement  it  is  necessary 
to  first  consider  the  system  adjoint  to  the  homogeneous  equations  (3.3).  This  adjoint 
system  consists  of  the  functions  F(i/’)  and  (j(i/»)  which  are  the  solutions  of  the  coupled 
equations 


d^G 

dip'^ 


fOo  ,  iM  , 

“  T - -  T - - - T 


ip  ~ 

t\ili^\  d^G 


I,  (f)  +  =  0, 


-f 


iXip^ 

~k^ 


k^ 


dip^ 

fQo 


_ 

V  P 


+ 


dG 


k'^ 


J  dv 


iM 

k^  ^  k^ 


+  1  )  G  -  p-  —  0, 


(3.7n) 


(3.76) 


with  the  associated  conditions  that  F,  G  and  dG/dtp  — >  0  as  ;/>  — >  0,oc.  Multiplying 
(3.6a)  by  F,  (3.66)  by  G,  adding  and  integrating  with  respect  to  xp  demonstrates  that 
for  a  solution  of  (3.6)  to  exist  necessitates  that 


where 


and  /*  =  Hdip,  where 


(3.Su) 

(3.S6) 


dU, 


20 


dp 


U^O  dV20 

2  dip 


-f  2F20 


dp'^ 


-  UMo  -  UooLho 


+ 


FFio^o 


G 

k 


dV^o 

dp 


+  Iv'o 


dV20 

dp 


dp  dp'^  2  dp  dp'^  dp^ 


V:,d^\o\ 
2  dp^  J 
(3.Sc) 
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Before  presenting  the  results  of  the  computations  we  observe  that  the  above  anal¬ 
ysis  may  be  easily  generalised  so  that  the  vortex  evolution  on  a  slow  lengthscale  and 
thence  the  growth  or  decay  of  a  near  neutral  disturbance  may  be  monitored.  For¬ 
mally,  if  we  consider  the  neighbourhood  of  the  point  X  —  Xn  (the  point  at  which  an 
infinitesimal  vortex  is  neutrally  stable)  and  allow  a  constant  frequency  disturbance  to 
develop  on  a  lengthscale  X  where  X  =  G~^ h~^(X  —  Xn)  with  h  1,  then  the  result 
of  repeating  the  previous  analysis  for  a  disturbance  of  amplitude  hA{X)  relative  to 
the  scalings  implied  in  (2.5)  is  the  evolution  equation 

It  is  then  simple  to  derive  the  equation  for  the  vortex  amplitude 

4-(|.4n  =cia,|.4|^-fC2|.4r,  (3.96) 

dX 

where 

Cl  =  2Re  ^  \  ~J  (3.9c,  d) 

In  order  to  evaluate  the  coefficients  in  (3.96)  it  is  necessary  to  solve  the  systems  (3.3-5) 
and  (3.7)  so  that  the  expressions  for  zi,  Z2  and  /,  as  defined  in  (3.8)  may  be  found. 
Several  limits  are  of  interest:  in  particular  we  discuss  asymptotic  explanations  of  the 
properties  of  solutions  of  equation  (3.96)  for  situations  of  large  or  small  scaled  vortex 
wavenumber  k  as  defined  in  (2.46)  and  (2.15a). 

It  is  clear  that  if  is  an  eigenset  of  (3.3)  then  so  is  (— —A)  for 

real  crossflows  A.  Evidently,  it  is  possible  to  restrict  our  attention  to  positive  crossflow 
parameters  without  any  loss  of  generality  and  this  is  done  for  the  remainder  of  the 
paper. 

§4  Numerical  methods  and  preliminary  calculations 

In  the  context  of  their  work  on  the  viscous  vortices  BH  were  solely  concerned  with 
the  numerical  solution  of  the  linear  system  (3.3).  Their  computations  were  based  on  a 
technique  originally  proposed  by  Malik,  Chuang  and  Hussaini  (1982)  in  which  the  dif¬ 
ferential  equations  to  be  solved  are  reduced  to  a  set  of  linear  algebraic  equations  using 
either  a  finite  difference  discretisation  or  a  spectral  representation  and  the  eigenvalues 
found  by  solving  the  characteristic  determinant  of  a  generalised  eigenvalue  problem. 
The  particular  method  described  in  Malik  et  al.  uses  a  fourth  order  accurate  (Euler- 
Maclaurin)  finite  difference  scheme  with  nodal  points  distributed  so  as  to  resolve  any 
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singular  layers  and  was  implemented  in  order  to  examine  the  temporal  and  spatial 
stability  of  a  three-  dimensional  compressible  boundary  layer  flow  over  a  swept  wing. 

For  the  present  work  we  attempted  to  adapt  the  methods  used  by  BH  but  even¬ 
tually  concluded  that  we  needed  to  use  a  more  efficient  algorithm  here.  The  rea;,ons 
for  reaching  this  conclusion  were  twofold:-  first  the  form  of  the  code  relevant  to  the 
linearised  problem  (3.3)  does  not  lend  itself  to  easy  modification  for  use  in  inlioinuge- 
neous  problems  such  as  (3.5).  Second,  and  perhaps  more  importantly,  BH  conducted 
a  small  number  of  computations  of  solutions  to  (3.3)  restricted  to  the  three  scaled 
frequency  choices  i7o  =  1,0, —1.  (This  was  due  to  limits  on  the  available  computing 
resources.)  In  order  to  execute  the  additional  calculations  required  to  compute  the 
amplitude  equation  coefficients  cq  and  C2  as  defined  in  (3.9)  and  to  extend  the  re:^'.il!r> 
to  a  greater  frequency  range  it  was  felt  that  an  improved  algorithm  was  required.  \\i- 
mention  in  passing  that  there  was  an  unexpected  by-product  of  using  the  improved 
scheme  as  it  enabled  us  to  explain  certain  of  the  computational  difficulties  encountered 
by  both  BH  and  Bassom  (1992)  and  which  are  described  later. 

We  first  detail  our  computational  techniciue  as  far  as  it  ;q:)plies  to  solving  the  ho¬ 
mogeneous  system  (3.3).  The  method  has  to  be  adapted  for  use  on  non-homogeneous 
problems  (for  example  (3.5))  and  these  modifications  will  be  described  as  reciuired. 
Therefore  we  note  first  that  (3.3)  can  be  characterised  by  the  pair  of  equations 

d^V  ,  dH'  .  dV  ^  , 

dy*  (lyi  —  0^  ld.l((j 

d-  d(uU  -h)5uF  =  0,  (d.li') 

where  the  coefficients  du,  d,,, . . .  etc.  are  known  functions  of  y. 

The  system  comprising  of  (4.1)  with  associated  boundary  conditions  typified  by 
(3.3c)  is  an  eigenvalue  problem  and  was  solved  by  considering  (4.1)  with  only  five  of 
the  homogeneous  boundary  conditions  invoked.  In  addition  a  normalisation  constraint 
was  imposed  so  that  (4.1)  was  solved  subject  to 

dU 

u  =  =  0,  —  I  at  y  =  0,  (.4.2u) 

dy 

U  =  V  =  - >  0  as  y  — >  oc.  (4.26) 

dy 

Iteration  on  the  eigenvalue  ensured  that  the  sixth  boundary  condition  dVjdy  =  U 
at  y  =  0  was  satisfied.  Specifically,  in  the  context  of  problem  (3.3),  system  (4.1  2) 
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was  solved  for  complex- valued  for  chosen  frequency  Qoi  crossflow  A  and  vortex 
wavenumber  k.  Iteration  on  the  imaginary  part  of  do  using  a  Newton  iteration  to  vary 
A  enabled  neutrally  stable  solutions  to  be  found:  cubic  splines  were  used  to  track  along 
the  neutral  curves. 

As  previously  mentioned,  the  main  benefit  of  the  present  algorithm  above  tliat 
used  by  BH  is  the  much  greater  speed  of  calculation  and  this  is  due  in  part  to  the 
implementation  of  a  variable  mesh  step.  The  equations  (4.1)  were  discretised  using 
central  differencing  on  the  grid  i/„,  n  —  1,2,  ...A’’,  where  2/1=0  and  yx  =  y^^K 
where  y^'^^  was  chosen  to  be  sufficiently  large  so  that  the  subsequent  results  were 
independent  of  its  value.  The  algorithm  permits  any  grid  discretisation,  although  it 
was  found  to  be  adequate  to  limit  the  form  to  four  distinct  regions  in  which  the  step 
length  was  constant:  specifically,  step  length  li  =  for  0  <  y  <  h  =  for 
<  y  <  y^-^K  j  =  2,3,4,  where  y^'^^  —  y^°^K  This  grid  was  found  to  both  resolve 
both  the  details  of  the  solution  and  also  integrate  over  the  entire  field  without  using  un¬ 
necessary  steps  in  the  outer  region.  In  all  the  calculations  described  henceforth  checks 
were  implemented  to  ensure  that  the  solutions  responded  smoothly  near  discontinuity 
points  in  the  mesh  -step  size. 

Typically  this  method  allowc'd  us  to  use  roughly  one-third  of  the  number  of  points 
as  was  required  by  BH  in  their  implementation  of  the  code  described  by  Malik  et  al. 
(1982).  Further  time  saving  was  obtained  by  formulating  our  code  so  as  to  account 
for  the  exact  structure  of  the  discretised  equations  to  be  solved.  The  method  used  in 
BH  may  be  applied  to  a  much  wider  variety  of  boundary  value  problems  whereas  the 
method  used  here  is  specific:  however  we  believe  that  restricting  ourselves  to  sucli  a 
system  specific  procedure  was  worthwhile  for  with  modest  grid  adaptation  we  achieved 
a  reduction  in  computer  time  to  roughly  one-thirtieth  of  that  needed  by  BH.  Details 
of  the  discretisation  used  and  of  the  procedure  for  solving  the  resulting  equations  are 
given  in  Appendix  A. 

As  a  check  on  the  method  described  above  it  was  decided  to  repeat  some  of  the 
linear  calculations  of  BH  to  ensure  consistency.  Although  BH  did  conduct  a  few  com¬ 
putations  relevant  to  non  neutral  infinitesimal  modes  it  should  be  remembered  that 
our  current  concern  is  ultimately  with  describing  the  weakly  nonlinear  modes  outlined 
in  section  3:  disturbances  which  by  definition  are  ‘close’  to  the  linear  neutral  stability 
curve.  Therefore  ('ffort  was  concentrated  on  tin.-  neutral  modes  of  BH  and  figure  1 
illustrates  such  modes  for  a  variety  of  scaled  frequencies  Dy.  BH  restricted  themselves 
to  the  choices  i'lo  =  1,0  and  —1,  but  the  efficienc}'  of  the  imjjioved  algorithm  allowi'd 
neutral  modes  to  be  derived  for  Do  =  —2(0. 5)3.  Figure  1  shows  the  crossflow  i)araiiu  u  r 
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A  and  figure  2  the  parameter  i3q  as  functions  of  the  scaled  vortex  wavenumber  A  .  W'e 
remark  that  our  calculations  faithfully  reproduced  the  corresponding  findings  of  Bli. 

Certain  trends  surmised  by  BH  on  the  basis  of  their  numerical  work  restricted  to 
the  three  values  of  Qq  are  confirmed  by  our  more  extensive  results.  In  particular  it  is 
seen  that  neutral  modes  appear  to  be  possible  over  wide  ranges  of  the  wavenumber  and 
that  for  large  k  the  crossflow  needed  to  produce  neutral  modes  is  quite  small.  Moreover. 
DHS  illustrated  that  for  flo  =  A  =  0  all  modes  within  the  0{G^  )  wavenumber  regime 
are  unstable;  consequently  crossflow  is  seen  to  have  a  stabilising  effect.  A  striking 
difference  also  occurs  depending  upon  the  sign  of  the  scaled  frequency  flu-  U 

the  neutral  modes  persist  over  th('  complete  range  of  A  and  as  A  — >  0  the  crossHow  A 
required  in  order  to  preserve  neutral  modes  tends  to  an  0(1)  value.  However  for  lly  <  b 
BH  and  Bassom  (1992)  could  not  find  such  disturbances  for  all  wavenumbers  A  and  the 
neutral  modes  apparently  disapp('arcd  at  some  critical  A.  Figure  1  demonstrates  that 
this  critical  value  increases  as  Qo  becomes  more  negative,  and  a  possible  explanation 
for  this  behaviour  is  given  later. 

For  small  wavenumbers  A,  figure  2  illustrates  that  the  behaviour  of  the  stream\'.-i.-?e 
variation  parameter  do  is  ^ilso  critically  dependent  upon  fly-  For  fly  >  0.  as  A  — -  0 
then  do  0  smoothly  (and  is  pro])ortional  to  A^  by  the  results  of  BH).  A  contrasting 
situation  occurs  for  flo  <  0  for  now  both  A  and  Ju  develop  erratic  belniviours  as 
A  decreases;  behaviour  which  was  also  found  in  BH.  Extensive  numerical  checks  on 
both  our  present  code  and  that  developed  for  use  in  BH  were  made  and  confirm  these 
findings. 

To  conclude  this  section  it  is  worthwhile  to  note  that  the  results  shown  in  figure 
1  are  sufficient  to  demonstrate  that  increasing  crossflow  stabilises  the  non-stationary 
vortex  modes.  There  are  several  regimes  in  which  asymptotic  description  is  possible 
(in  particular  those  cases  of  large  and  small  A)  but  since  the  primary  aim  of  the 
present  study  is  to  examine  the  weakly  nonlinear  system  we  defer  these  asymptotic 
considerations  for  the  present. 

§4.1  The  weakly  nonlinear  calculations 

Given  the  outline  of  the  nuiiK'rical  method  above  the  imi)leinentation  of  the  rou¬ 
tines  required  to  evaluate  the  aini)litude  equation  coefficients  ci  and  cj  was  straightfor¬ 
ward  and  the  following  methodology  was  adopted.  Once  the  homogeneous  ecpiatioiis 
(3.3)  were  solved  for  the  streamwise  variation  parameter  ,du  and  scaled  crossflow  A  in 
terms  of  the  given  frequency  flo  <’iid  wavenumber  A  and  the  resjx'ctive  eigenfunctions 
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L'^io  and  V'lo  deduced  it  was  a  simple  matter  to  evaluate  the  mean  flow  quantities  Uuu 
and  $0  deflned  in  (3.4). 

Next  the  solution  of  the  inhomogeneous  system  (3.5)  was  required  in  order  to 
determine  the  second  harmonic  terms  U20  and  Appendix  A  details  the  solu¬ 

tion  strategy  for  the  homogeneous  problem  but  certain  modiflcations  are  required  for 
dealing  with  an  inhomogeneous  case.  These  tire  briefly  described  in  Appendix  B.  Fi¬ 
nally,  the  adjoint  functions  F(i^’)  and  G(</’)  as  given  by  (3.7)  were  found  by  using  the 
ideas  presented  in  Appendix  A.  The  integrals  ci  and  C2  as  deflned  in  (3.8)  were  then 
evaluated  using  Simpson’s  rule  and  the  amplitude  equation  coefficients  Cj,  c>  thence 
deduced.  The  usual  checks  were  made  to  ensure  the  accuracy  of  the  numerical  solutions 
of  the  systems  (3.3),  (3.5)  and  (3.7).  In  addition,  one  further  check  was  implemented 
for  each  set  of  solutions;  in  each  system  the  discretised  solutions  were  substituted  into 
each  side  of  the  equations  to  cn.'^ure  that  the  required  matrix  equations  were  being 
solved  accurately  and  that  the  derived  results  did  indeed  solve  the  problem. 

We  now  consider  some  properties  of  the  resulting  amplitude  equation  for  a  variety 
of  parameter  choices. 

§5  Remarks  on  the  weakly  nonlinear  properties  of  the  vortex  modes 

Implementation  of  the  numerical  procedures  described  in  section  4  led  to  the 
determination  of  the  coefficients  C]  and  C2  appearing  in  the  amplitude  equation  (3.96). 
This  equation  demonstrates  that  the  weak  nonlinearity  of  the  problem  allows  the 
existence  of  a  threshold  equilibrium  amplitude  .4c  given  by 

|.4cP  =  -^fh2,  (5.1) 

with  fl2  >  0  or  1^2  <  0  as  appropriate  in  order  to  ensure  that  the  right  haiul  sitle  of 
(5.1)  is  positive. 

Calculations  were  conducted  for  a  variety  of  scaletl  frequencies  Qy  taking  values 
between  —2  and  3  and  the  results  are  sunnmirised  in  figures  3-5.  Figure  3  illustrates 
the  dependence  of  the  coefficient  of  the  linear  term  in  (3.96)  (i.e.  Cj)  upon  flo  and 
the  vortex  wavenumber  k.  It  is  convenient  to  discuss  the  cases  Qu  ^  0  and  Qy  <  0 
separtitely.  In  the  former  case  ci  >  0  across  the  whole  range  of  wavenumber  sptice 
which  indicates  that  according  t(}  linear  theory  a  mode  of  frequency  Qu  +  is 

unstable  for  Q2  >  0  whereas  it  is  stable  for  Q2  <  0.  This  is  to  be  expected  as  from 
figure  1  it  may  he  deduced  that  for  a  fixed  crossflow  then  as  the  frequency  of  the  mode 
increases  so  it  loses  stability.  Furthermore,  figure  3  also  demonstrates  that  for  flo  <  U 
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a  more  erratic  behaviour  for  ci  occurs.  Specifically  there  are  regions  of  wavenumber 
space  where  ci  just  becomes  negative  and  there  are  rapid  variations  in  ci  for  small 
k.  These  changes  are  easily  explained  by  an  asymptotic  analysis  of  the  governing 
systems  (3.3),  (3.5)  and  (3.7)  for  small  k:  an  analysis  which  is  closely  allied  to  work 
elucidated  to  in  BH  and  which  is  explained  in  section  G.  However  the  small  k  case  for 
fio  <  0  evolves  to  a  very  different  structure.  Although  BH  postulated  upon  the  form 
of  this  revised  configuration  they  could  not  solve  the  resulting  equations-  the  reason 
for  this  was  surmised  by  Bassom  (1992)  as  being  due  to  a  cut-off  point  beyond  wliicli 
no  neutral  solutions  could  exist.  The  curves  of  figure  3  for  flu  <  0  lend  some  credence 
to  this  suggestion  for  as  Qq  becomes  increa.singly  more  negative  so  the  calculations 
become  increasingly  more  difficult  to  perform  for  snudl  k.  However,  we  have  also  some 
evidence  that  a  much  more  elegant  computational  procedure  may  lead  to  a  resolution 
of  these  difficulties;  a  matter  which  is  discussed  at  more  length  later. 

A  final  comment  on  figure  3  concerns  the  wild  ‘dips'  in  cq  observed  at  relativel}^ 
large  values  of  k  when  CIq  <  0.  This  feature  is  not  seen  for  positive  Qy  and  extensive 
numerical  checks  of  our  work  have  suggested  that  this  behaviour  is  real  and  not  just 
some  numerical  artefact.  As  yet  we  have  no  convincing  argument  as  to  why  this 
behaviour  should  occur. 

We  turn  now  to  interpret  the  findings  summarised  in  figure  4  whicli  illustrates 
the  dependence  of  coefficient  C2  on  k.  This  coefficient  is  of  greater  importance  than 
Cl  for  the  sign  of  C2  reveals  the  nature  of  the  equilibrium  amplitude  as  defined  in 
(5.1).  For  all  frequencies  fly  >  0  it  is  seen  that  c-z  <  0  and  consequently  the  weak 
nonlinearity  of  the  flow  leads  to  a  stabilisation  of  the  vortices.  In  this  case  vortices 
with  initial  amplitude  less  than  Ae  grow  to  Ae  whereas  those  with  initial  size  gretiter 
than  Ae  tend  to  diminish.  However  it  must  be  recalled  that  if  the  scaled  amplitude  A 
becomes  too  large  the  fundamental  assumptions  underlying  weakly  nonlinear  theory 
are  invalidated  and  a  fully  nonlinear  account  of  the  vortex  structure  is  essential. 

In  the  main  the  finding  that  C2  <  0  also  holds  for  negative  scaled  frequencies 
The  one  additional  feature  that  appears  now  is  the  occurence  of  a  large  negative  spike 
in  the  value  of  C2  which  is  evident  at  a  relatively  small  value  of  k.  .4gain  extensive 
numerical  checks  suggest  that  this  is  indeed  a  true  feature  as  opposed  to  a  computa¬ 
tional  peculiarity.  Figures  4e-g  also  demonstrate  that  the  size  of  this  spike  decreases 
as  fly  becomes  increasingly  more  negative. 

Finally  for  this  section,  equation  (5.1)  trivially  yields  information  concerning  the 
equilibrium  amplitude  Ae  and  which  is  shown  in  figure  5  for  the  various  values  of 
fly.  Of  course  many  features  of  this  amplitude  follow  directly  from  the  information 
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portrayed  in  figures  3  and  4  and  so  require  little  additional  comment.  We  mention 
first  that  for  all  the  cases  considered  C2  <  0  so  that  modes  with  perturbed  frequency 
^2  >  0  would  be  expected  to  evolve  to  these  equilibrium  amplitudes  whereas  those 
with  fl2  <  0  would  decay  to  zero.  Again  it  is  convenient  to  divide  our  discussion 
between  the  cases  Qq  >  0  and  fio  <  0-  In  the  former,  as  A:  — >  oo,  j-4e)  ^  0  and  as  k 
decreases  \Ae  \  rises  and  reaches  a  maximum  value.  As  Qq  increases  so  this  maximum 
value  decreases.  The  form  of  figures  5a-c  suggest  that  for  vortices  of  scaled  frequency 
flo  >  0  those  modes  with  low  wavenumbers  relative  to  the  O(G^)  implied  scaling 
would  be  more  readily  observed  than  those  of  greater  wavenumber.  It  should  also  be 
noted  that  as  A:  0  the  calculations  became  increasingly  sensitive  to  the  grid  spacing 
used  and  the  slight  wobbles  in  the  curves  as  A*  — >  0  are  a  direst  consequence  of  this. 

For  Qq  negative  a  number  of  the  trends  just  described  continue  to  persist.  In  pai- 
ticular,  as  fio  moves  through  increasingly  more  negative  values  |Ae|  decreases  although 
its  maximum  value  for  any  fixed  frequency  occurs  at  increasing  values  of  k.  Again  dif¬ 
ficulties  in  evaluation  were  experienced  for  small  wavenumbers  and  furthermore,  the 
small  values  of  C2  as  demonstrated  in  figures  4e-g  make  computation  of  |.4e|  at  large 
k  difficult. 

To  conclude  this  description  of  the  weakly  nonlinear  calculations  it  is  worthwhile  to 
briefly  recall  the  more  salient  results.  Importantly,  whereas  BH  showed  that  the  effect 
of  crossflow  is  to  stabilise  the  OfGs)  wavenumber  modes  according  to  a  linearised 
theory,  the  above  work  has  demonstrated  that  crossflow  also  stabilises  the  modes 
on  a  weakly  nonlinear  basis.  The  corresponding  equilibrium  amplitudes  calculated 
over  a  range  of  frequencies  and  wavenumbers  reveals  that  the  largest  amplitudes  are 
associated  with  near  stationary  vortices  and  it  can  be  tentatively  proposed  that  these 
modes  are  the  most  likely  candidates  for  practical  observation.  The  analysis  origiirally 
proposed  by  BH  for  small  wavenumber  linearised  vortices  at  scaled  frequencies  f7o  >  0 
may  be  adapted  to  allow  discussion  of  these  modes  using  a  weakly  nonlinear  approach 
and  this  is  considered  now. 

§6  The  low  wavenumber  (A-  <C  1)  limit  for  flo  >  0 

The  calculations  presented  in  BH  showed  that  when  Qq  >  0,  A’  <C  1  the  eigenfunc¬ 
tions  Uio  and  Fio  of  system  (3.3)  are  concentrated  in  a  thin  region  near  0  =  0.  They 
demonstrated  that  these  functions  assume  a  multi-zoned  structure  form  and  that  the 
crossflow  A  and  streamwise  variation  parameter  /^o  required  for  neutral  modes  assume 
the  forms 

A  =  Ao  -f  A]  A’^  ,  Po  —  /^ooA^  -t-  /^oi  A^  -f-  . . . ,  (G-l) 
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where  all  the  constants  are  real  valued.  The  asymptotic  structure  for  the  solution  is 
summarised  in  figure  6  where  it  is  shown  that  the  configuration  divides  into  a  main 
zone  {II  on  figure  6)  of  thickness  0{k^ )  which  contains  a  thin  region  of  depth  0{k^ ), 
I,  supplemented  by  a  viscous  wall  layer  of  thickness  0(A.-),  Ill,  and  a  far  field  zone, 
IV.  It  is  to  be  recalled  that  the  normalisation  chosen  for  all  the  numerical  work  in 
this  paper  is  that  U[q  =  1  at  ^  =  0.  Guided  by  this  requirement  it  is  easy  to  adapt 
the  workings  of  BH  to  show  that  the  most  important  contributions  to  the  integrals 
defined  in  equations  (3.8)  occur  within  the  zone  I  wherein  it  is  convenient  to  write 

t/)  =  k'^il'o  +  (^const.  +  ifj  .  (C.2a) 

3 

Here  ipo  ~  1.47fio  and  the  constant  has  a  value  which  is  of  no  concern  for  the  current 
workings.  Furthermore,  in  zone  I  the  solutions  Uiq  and  Fio  develop  according  to 

U\Q  =  k  4^10(1^)  +  ...,  Fio  =  Fio(^)  +  •  ■  • )  (G.2i) 

where  these  functions  satisfy 

(  <P  .  ' 

I  ~  tAoV* 

\dip 


(6.2d) 


(6.2e) 

The  constant  C  is  chosen  so  that  U[q  =  1  at  0  =  0  and  elementary  matching  between 
the  solutions  valid  in  each  of  the  regions  sketched  in  figure  6  leads  to 

C«-0.21iQj.  (6.3) 

The  eigenproblem  (6.2)  was  first  solved  by  Hall  (1985)  who  showed  that 

0oAo  ^  w  4.69,  ^  -2.89, 
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which  in  turn  gives  Aq  ^  0.46fio  .  This  prediction  of  crossflow  required  in  order  to 
produce  neutral  modes  as  ^  >  0  gives  good  agreement  with  the  forms  depicted  in 

figure  1  for  frequencies  fio  >  0. 

Following  this  work  taken  from  BH  it  is  straightforward  to  deduce  the  forms  of 
the  remaining  functions  within  the  thin  zone  I.  The  adjoint  functions  are 


where 


F  =  F(V0  +  ...,  G  =  kG{^)  +  ..., 


d?  -  c. 

-  iQ.  -  iXo-ip  I  F  +  i>oG  -  0, 

dxp  J 


(P  .A  }dG  ^ 

—  til  —  iXqxJ)  I  — ^  =  r  , 

dip  j  <ii/> 


subject  to  the  conditions  that 


(G.4a) 


(6.46) 


F  cx  —  and  G  oc  as  j-^l  — >  oo. 

t/>  rp 


Furthermore,  U2Q  =  k  *U2o  +  •••,  V20  =  ^'•*^20  +  •••  where  U20  and  F20  are  odd¬ 
valued  functions  satisfying 


—  2zS2  —  2iAQ'tl>  j  1/20  —  V20  —  — ^1^10  “  U 10 — — 
dxp  j  dtp  d^ 


(6.5a) 


d^  -  "^\d^V  ~  I 

— — ‘liCl  —  2iXQ'ip  j  — ^2 — h  4fAoF2o  4"  4^0^20 

dip  J  dip 

=  2 


^dViod^Vu 

0,3  “  "X  -2 

dip  dip  dip 


(6.56) 


subject  to  the  constraints  that  U20  =  0{ip  )  and  F20  =  0(ip  )  as  |0|  — >  cx).  Finally, 
within  zone  I  the  functions  Uoo  and  4»o  are 


b’oo  —  k  *U 00(0)  +  •  •  • ,  ^0  =  ^oo(0)  +  •  •  • , 
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where 


U 


00 


■/:.( 


U loV  10  +  U iqV  10  )  d4^  and  $ 


■00 


-i: 


V 


dV 


10  ■ 


10 


d4) 


d4’. 


(6.G) 

Formally,  it  is  possible  to  use  the  expansions  given  above  to  deduce  the  forms  of  the 
various  flow  quantities  in  each  of  the  remaining  regions  II-IV .  However  tedious  but 
straightforward  manipulations  verify  that  the  dominant  contributions  to  the  expres¬ 
sions  (3.8a  —  c)  arise  from  within  zone  I. 

The  systems  (6.2),  (6.4)  and  (6.5)  were  solved  using  techniques  based  upon  the 
general  descriptions  given  in  Appendices  A  and  B  (some  additional  remarks  specific 
to  these  systems  are  also  in  the  latter  Appendix).  Simpson’s  rule  was  used  to  evaluate 
the  necessary  integrals  and  thence  the  amplitude  equation  coefficients  were  found  to 
take  the  forms 


C\  %  0.7f7g  ®  -f- . . . ,  C2  ~  — 7^0**  -f- . . . ,  (6. 1  Cl.  b) 

as  k  0.  These  asymptotic  predictions  provided  reasonable  agreement  with  the 
numerical  results  discussed  in  the  preceding  section.  Results  (6.7)  then  imply  that  the 
equilibrium  amplitude 

|Ae|  «  0.3^0“ 

which  again  suggests  that  it  is  the  low  wavenumber  relatively  low  frequency  modes 
which  have  the  greatest  equilibrium  amplitudes.  However  there  is  evidence  from  figures 
5a-c  that  as  k  decreases  the  accuracy  of  our  computational  work  deteriorates.  One 
difficulty  with  the  prediction  of  Ag  as  k  0  is  provided  by  the  extreme  smallness 
of  the  coefficient  C2  (see  (6.76)).  Clearly  tiny  inaccuracies  in  the  evaluation  of  co 
can  have  drastic  consequences  for  the  resulting  value  of  |.4e|.  As  A-  — >  0  it  has  been 
demonstrated  that  the  entire  solution  structure  becomes  compressed  against  the  wall 
and  then  the  numerical  resolution  of  the  distinct  zones  sketched  in  figure  6  is  rendered 
very  difficult.  Limits  on  our  numerical  work  were  imposed  by  the  available  resources 
and  it  was  considered  that  it  would  not  be  worthwhile  to  attempt  to  refine  our  grid 
further  for  reasons  described  presently. 

The  scalings  described  above  fail  as  Qg  ^  0  foi’  this  limit  the  thin  layer  / 
moves  towards  the  wall  and  Bassom  (1992)  has  demonstrated  that  when  flo  =  0{k^) 
this  layer  merges  with  the  viscous  wall  layer  III.  Within  the  wall  layer  the  governing 
equations  now  take  the  forms 


(i  -  ,!i  -  -  a?)  iJ  =  V. 


{G.Sci) 
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(6.86) 


f  -  zfi  -  iM  -  ixiA  ^  +  2i\v  =  -ii), 

J  (le 

where  the  *  denotes  quantities  scaled  on  suitable  powers  of  k  (fuller  details  are  given 
in  Bassom  (1992)).  In  particular,  Qq  =  and  these  wall  layer  equations  need  to  be 
solved  subject  to  the  conditions  that 

~  ~  dV 

U,  V  and  — ^  — *  0  at  ^  =  0  and  as  ^  oo.  (6.8c) 

d^ 

Bassom  (1992)  used  the  numerical  code  of  Malik  ei.  al.  (1982)  to  solve  (6.8) 
for  a  selection  of  values  of  Q.  As  11  — >  oo  he  showed  how  these  low  frequency  modes 
match  with  those  for  which  Ho  =0(1)  but  he  could  only  manage  to  obtain  solutions 
of  (6.8)  for  H  greater  than  He  ~  —3.3.  The  work  of  Bassom  (1992)  was  prompted  in 
part  by  the  finding  of  BH  concerning  the  low  wavenumber  structure  outlined  by  (6.2) 
above.  The  asymptotic  regions  found  bj'^  BH  are  only  valid  for  positive  v^alues  of  Ho 
and  that  paper  made  tentative  suggestions  concerning  possible  configurations  when 
Ho  <  0.  In  particular,  a  critical  layer  type  problem  was  proposed  but  no  attempt 
was  made  to  solve  this.  Subsequent  calculations  failed  to  find  a  solution  and  Bassom 
(1992)  hoped  to  connect  the  Ho  >  0  and  Ho  <  0  cases  by  investigating  the  intermediate 
problem  where  Ho  =  k'^H  by  examining  the  numerical  solutions  of  (6.8)  as  H  — >  —  oc. 
However  the  existence  of  the  cut-off  frequency  He  below  which  numerical  solutions 
were  unattainable  thwarted  this  aim. 

The  system  (6.8)  was  solved  using  slight  modifications  to  the  numerical  techniques 
described  in  Appendix  A.  Like  Bassom  (1992),  difficulties  were  encountered  as  He 
was  approached.  Closer  studies  of  this  phenomenon  revealed  that  the  as  H  —>  He 
the  eigensolution  migrates  from  the  wall  ^  =  0  and  thus  becomes  almost  completely 
independent  of  the  boundary  conditions  imposed  at  the  wall.  Therefore  the  solution 
becomes  insensitive  to  the  values  of  the  eigenparameters  and  numerical  convergence  is 
impossible  to  obtain.  This  type  of  behaviour  typically  occurs  whenever  a  ‘null  space’ 
of  the  system  is  approached  and  indicates  that  the  elimination  techniques  employed 
to  solve  the  discretised  equations  need  to  be  replaced  by  a  scheme  which  directly 
inverts  the  entire  discretised  set  (such  schemes  are  termed  ‘global  methods’).  This  is 
computationally  extremely  expensive  but  it  does  lead  to  a  determination  of  a  spectrum 
of  eigenvalues  of  the  problem. 

This  difficulty  which  arises  when  a  null  space  of  a  problem  is  encountered  also 
provides  the  explanation  for  some  of  the  other  deficiencies  in  our  numerical  work  to 
date.  In  figures  1  and  2  where  neutral  curves  were  presented  it  is  noted  that  as  Ho 
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becomes  progressively  more  negative  the  wavenumber  range  over  which  results  were 
obtained  diminishes.  Investigations  have  shown  that  this  is  again  due  to  approaching 
a  null  space  of  the  appropriate  system  (3.3). 

BH  also  made  some  remarks  concerning  the  properties  of  large  wavenumber  modes. 
They  demonstrated  that  as  k  oo  the  disturbance  remains  confined  within  an  0(1) 
thick  region  which  moves  away  from  0  =  0.  Specifically,  if  0  =  ck^  +  const.  +  0^ 
(where  the  precise  value  of  the  constant  is  of  no  consequence  for  the  leading  order 

problem),  /3o  =  l^ook^  +  . . .  and  A  =  Xok~^  + _  If  the  eigenfunctions  of  (3.3)  are 

written  as 

Uio  =  u}  +  ...,  V^o  =  kX^+...,  (6.9a) 

then  it  is  easily  verified  that  ^oo  =  —  Aqc  and  the  eigenfunctions  satisfy 


T/t  _  _rTJ^ 

Kq  —  CUq  , 


(6.96) 


(  ^  -  1  +  j  0  =  C.t .  (6.9c) 

Again  BH  were  unable  to  solve  this  eigenproblem  for  ^oo  and  c  and  it  is  observed 
from  that  paper  that  calculations  of  the  homogeneous  system  (3.3)  failed  for  values 
of  k  larger  than  about  1.4.  Further  researches  have  once  more  confirmed  thai  these 
problems  are  due  to  approaching  a  null  space  of  equations  (3.3). 

In  order  to  develop  our  findings  of  this  paper  further  work  is  needed  to  investigate 
the  global  numerical  methods.  Preliminary  runs  with  a  small  number  of  grid  points 
have  been  conducted  for  the  eigenproblem  (6.2)  first  solved  by  Hall  (1985)  and  for 
the  small  frequency  problem  (6.8)  first  written  down  by  Bassom  (1992).  Results  for 
this  latter  problem  have  revealed  the  existence  of  solutions  for  Cl  less  than  the  cut-off 
value  Clc.  This  indeed  proves  that  the  previous  problems  for  <  fic  are  entirely 
numerically  based  and  have  no  physical  significance.  As  yet  there  are  too  few  results 
from  the  computationally  intensive  global  method  in  order  to  justify  a  fuller  discussion 
here.  However  it  probable  that  developments  with  these  ideas  will  lead  to  an  improved 
explanation  of  the  large  wavenumber  and  negative  flo  limits  and  we  plan  to  report  on 
our  findings  in  due  course. 


§7  Discussion  and  Conclusion 

In  the  previous  sections  we  have  described  the  mechanism  by  which  the  intro¬ 
duction  of  crossflow  influences  the  weakly  nonlinear  stability  characteristics  of  0(Gs  ) 
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wavenumber  viscous  vortices.  The  work  was  confined  to  the  situation  in  which  the 

_  Jl. 

crossflow  is  of  size  0(i?e  ^ )  since  that  is  the  crucial  size  at  which  the  crossflow  first  lias 
a  significant  effect  upon  on  the  vortices.  Additionally  the  Gdrtler  number  was  taken 
to  be  large  as  this  allows  non-parallel  effects  to  be  accounted  for  using  asj'inptotic 
means  and,  furthermore,  is  relevant  for  many  realistic  boundary  layer  flows. 

One  advantage  of  concentrating  on  the  0(G^  )  wavenumber  modes,  apart  from  the 
fact  that  these  are  the  most  unstable  for  a  two-dimensional  boundary  layer,  is  that  tlic 
structure  of  the  disturbance  permits  all  the  basic  flow  quantities  which  are  functions 
of  the  particular  boundary  layer  to  be  scaled  out  of  the  problem  leaving  a  system  of 
equations  which  is  valid  for  a  wide  variety  of  three-dimensional  flows.  Consequently 
the  scalings  (2.15)  would  need  to  be  applied  again  in  ortlcr  to  assess  the  implications 
of  our  findings  for  any  specified  flow. 

The  main  conclusion  of  our  work  has  been  that  over  all  the  wavenumbers  and 
frequencies  investigated  the  influence  of  weak  nonlinearity  is  stabilising  and  the  con¬ 
sequent  supercritical  equilibrium  amplitudes  have  been  evaluated.  These  amplitudes 
tend  to  be  largest  for  vortices  of  small  wavenumber  and  frequencies  relative  to  the 
initial  scalings  and  suggest  that  it  is  these  modes  which  would  appear  to  be  the  most 
likely  to  be  observed  in  practice.  A  limited  asymptotic  study  has  been  accomplished 
which  yields  indications  of  the  solution  characteristics  for  small  k  and  positive  scaled 
frequencies  flo-  However  the  corresponding  work  for  negative  Qq  ‘^nd  that  relevant  to 
large  wavenumbers  k  was  not  comjjleted  due  to  difficulties  in  encountering  null  spaces 
of  the  governing  differential  systems.  As  reported  in  the  previous  section  preliminary 
work  has  begun  using  more  appropriate  numerical  methods  in  order  to  circumviut 
these  problems.  We  feel  that  these  calculations  are  important  for  the  following  rea¬ 
son.  The  work  described  by  BH  and  that  here  has  conclusively  demonstrated  that 
both  the  linear  and  weakly  nonlinear  properties  of  these  viscous  vortices  are  critically 
dependent  upon  the  sign  of  the  scaled  frequency  and  crossflow  A.  For  any  particular 
boundary  layer  either  case  may  be  the  more  relevant  (depending  upon  the  signs  of  the 
scaling  quantities  within  (2.15))  and  so  it  is  desirable  that  both  eventualities  are  anal¬ 
ysed  properly.  Whereas  the  solution  properties  are  now  reasonably  well  understood 
for  Qo  >  0  this  understanding  is  clearly  deficient  for  Qo  <  0  ‘i^d  work  on  this  latter 
case  is  continuing. 

The  other  extension  of  this  paper  which  merits  close  study  is  the  question  as  to 
the  nature  of  the  fully  nonlinear  properties  of  viscous  vortices  in  three-dimeusiuiml 
boundary  layers.  The  investigation  of  Denier  &  Htdl  (1991)  described  in  the  introduc¬ 
tion  has  concentrated  upon  the  two  dimensioind  version  of  this  problem  and  slu)weil 
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that  in  general  the  vortex  suffers  a  finite  distance  algebraic  breakdown  as  it  develops 
downstream.  Denier  &  Hall  were  unable  to  conduct  a  weakly  nonlinear  analysis  of  the 
type  discussed  here  for  in  the  two-dimensional  boundary  layer  there  are  no  neutral 
modes  of  wavenumber  0(Gs).  The  results  of  BH  and  the  extensions  presented  in 
the  current  work  show  that  the  effect  of  three-dimensionality  is  to  stabilise  the  vortex 
mode.  Therefore  the  three-dimensional  version  of  Denier  Hall  (1991)  would  be  most 
valuable  in  deciding  whether  fully  nonlinear  or  crossflow  effects  are  the  more  impor¬ 
tant.  The  former  induces  a  catastrophic  breakdown  in  the  flow  whereas  the  latter  effect 
is  stabilising  and  so  it  can  be  anticipated  that  there  may  well  be  a  delicate  balance 
between  the  two.  It  is  noted  that  the  full  three-dimensional  nonlinear  equations  have 
already  been  obtained  (in  (2.16))  and  work  is  underway  to  investigate  the  properties 
of  their  solutions. 

In  the  context  of  our  work  concerning  the  viscous  modes  it  should  be  remembered 
that  a  concurrent  study  by  Blackaby  Sz  Dando  is  dedicated  to  examining  the  nonlinear 
properties  of  0(1)  wavenumber  inviscid  modes  whose  linear  stability  properties  were 
outlined  in  BH.  As  yet  we  are  not  aware  of  any  results  arising  from  their  work.  Their 
findings  will  need  to  be  compared  to  ours  in  order  for  the  relative  importance  of  the 
two  nonlinear  mode  types  to  be  accurately  assessed.  BH  attempted  to  resolve  the 
question  as  to  which  mode  is  the  more  likely  candidate  for  practical  observation.  B\- 
carrying  out  a  linearised  receptivity  calculation  of  the  type  given  in  DHS  relevant  to 
Gortler  vortices  in  two-dimensional  boundary  layers,  BH  were  able  to  show  that  wall 
roughness  is  a  more  efficient  stimulator  of  the  viscous  modes  than  the  inviscid  ones  and 
thus  the  viscous  modes  might  be  the  easier  to  generate  experimentally.  However  they 
also  pointed  out  that  as  the  crossflow  increases  the  growth  rates  of  the  inviscid  modes 
increase  to  become  larger  than  the  viscous  rates.  Thus  beyond  a  certain  crossflow  size 
the  observed  instability  may  well  be  a  Rayleigh  instability  as  opposed  to  a  centrifugal 
one. 

In  many  practical  situations  where  Gortler  vortices  are  thought  to  be  a  likely  cause 
for  transition  the  basic  state  is  three-dimensional.  Our  work  has  demonstrated  that 

_  i.  3  ^ 

a  crossflow  of  small  size  0{Re  is  sufficient  to  stabilise  vortex  modes  according 

to  a  weakly  nonlinear  basis  but  the  results  of  Denier  &  Hall  (1991)  indicate  that  full 
nonlinearity  of  the  disturbance  may  well  lead  to  rapid  breakdown.  It  is  a  matter  of 
some  interest  as  to  the  relative  imi)ortancc  of  these  crossflow  and  nonlinear  mechanisms 
which  can  only  be  resolved  by  extending  our  findings  to  the  fully  nonlinear  regime.  This 
problem  should  obviously  be  a  to])ic  of  careful  theoretical  and  practical  investigations 
of  the  Gortler  mechanism  in  thre<'  dimensional  boundary  layers. 
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APPENDIX  A 

In  this  appendix  we  briefly  describe  our  solution  strategy  for  equations  of  the 
general  form  (4.1). 

The  derivatives  in  (4.1)  may  be  discretised  as, 


_  ^n  +  l  4^n  —  \ 

dy  yn+l  -  Vn-l  ' 


(-41) 


(p<j) 

dy'^ 

and 


2(f>n 


+  1 


2(pn-\ 


{yn+i  -  y„-i)(y„-(-i  -y„)  (y„_i  -y„)(y„.^i  -  y„)  (y„+i  -y„_i)(y„  -y„_i)’ 

(-42) 


+2 


dV  _ 

dy^  (yn-2  -  y„-|-2)(yn-)-l  -  y„-|-2)(yn+2  -  yn-\){yn+2  -  yn) 

24<^„-i-i 


+ 

+ 

+ 


(yn-2  -  yn-|-l)(yn  +  l  -  yn-l)(yn  +  2  -  yn  +  l)(yn-|-]  -  yn) 

_ 24^n _ 

(yn— 1  ~  yn)  (yn— 2  yn)  (yn  yn+1 )  (yn  ~  yn+2) 

_ 24^n-l _ 

(yn-2  -  yn-l)(yn+l  “  yn-l)(yn+2  “  yn-l)(yn  ”  yn-l) 

_ 24^n-2 _ 

(yn-2  -  yn-l) (yn+1  “  yn-2)(yn+2  “  yn-2) (yn-2  “  yn)' 


(A3) 


Expressions  (Al-3)  were  all  obtained  by  implementing  the  symbolic  manipulation  pack¬ 
age  Mathematica  and  enable  the  homogeneous  equations  (4.1)  to  be  written  in  the 
general  forms 


^n ^n  +  2  4"  bfi Un+ 1  A  4"  d^l^n  —  1  4"  Cn ^n— 2  4"  fn —  4^n  >  (^  —  ‘  ‘  ^  I )  (-44ci) 

yn^n+1  4"  bfiUfi  -f-  ifiUfj  —  i  -j-  ~  Un-  (u  =  2,  ■  ■  ■  1)  (A46) 

For  the  homogeneous  problem  V),  =  =  0  for  all  n  and  the  boundary  conditions  are 

implemented  by  imposing 


Ui  =  t’l  =  0,  uyv  =  un  =  0) 

=  (Normalisation)  ""f  - =  Q.  (.«) 

y2  -  yi  y/v+i  -  yN 

In  the  context  of  system  (3.3)  it  remains  to  ensure  that  dVio/dy  =  0  at  the  wall  by 
varying  the  parameter  0.  The  algorithm  employed  to  solve  (.44,5)  was  an  essentially 
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standard  elimination  technique.  It  is  convenient  to  introduce  the  notation  Ev^j  and 
Ey,,j  to  refer  to  (A4a)  and  (A4b)  respectively  when  n  =  j . 

The  first  step  is  to  impose  the  lower  boundary  conditions  in  Eu,2  by  redefining 
U2  according  to 

U2  —*  U2  ~  ^2^2, 

where  U2  =  u\  +  {y2  —  y\)  u'  =  y2  —  yi,  since  a  =  0  and  a'  =  1  at  y  =  0.  Relations 
Eu,2  and  Ey^z  are  then  solved  for  V2  and  as  and  this  solution  is  written  as 

V2  =  (1v,2V2  +  ^v,2'^i  +  +  Rv,2i 

where 


^u,3 

with  Fs  = 

Uz  Uz  —  *3^2) 


U3  =  0.u,3Vz  +  iu,3*^4  +  <^u,3^5  +  -Ru,3) 


V.  -  EEx 

97 


=  ^3  f  ^3  -  ^3  ^  *^3 

2  ^  ^  ^i>,2  —  1^  ^  ^  7 

’  Ta  ’  Fa  Fa  la 

i2«v,2  f  J2bv,2  =  j‘2Cv,2  ^  U2  j2Ru,3 

—  ~  1  0u,3  —  7  Cti,3  —  ,  tlu,3  — 

92  92  92  92  92 

=  dz  —  f3j2/92-  Eu,3  is  considered  next  and  upon  redefining 


and  eliminating  uz  it  is  found  that 

93U4  +  j3Vz  +  n3a4  +  jdzVh  =  Uz, 
where  use  has  been  made  of  the  transformations 

j3  — >  is  +  hzOy^z,  oiz  =  hzbu^z,  ^3  —  hzCy^z,  Uz  — +  Uz  —  hzRu,3- 
Meanwhile  V2  can  also  be  eliminated  from  by  implementing 

^4  — ^  (^4  -f  ^40v,2,  C4  — >  C4  +  646,7,27  64  — >  64  +  646,7,27  V4  — >  V4  —  64i?,,,2- 

At  this  point  it  is  possible  to  write  this  procedure  in  terms  of  the  general  equations 
Ey^n  and  Solving  for  a„_2  and  u„_i  from  Ey^n-i  and  Ey^n-2,  gives 

Vfi  —  2  —  Ov,n  —  2^n  —  l  4"  6p,n  — 2^’m  4"  6y,n  — 2t’ri-(- 1  4"  Rv,n  —  2,  (.46fl) 

U|,  — 1  =  6tu,n  — l^n  — 1  4"  6u,n  — l^n  4*  6u,,j_ian4-l  4"  Ru  ,n  —  \  ,  (A66j 
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where 


r*  /n  —  1  Qm  —  a  J  /n  - 1  —  2 

=  9n-j  f  _ 

Ow,n— 2  —  ;;  )  Ot)^Tt—2 - 

1  n-1  1 n-1 

Vfn-lUn-i 

-:;  _  -I.  —  *  £)  _  ffn-2 

2  —  -p  )  ■^v,n—2  —  p  > 

■•■  n  — 1  1  n  — 1 


=  _  ^n  — 2  jn  —  20-v,n—2  t  ^n  —  2  jn  —  2^v,n—2 

Ou,n  — 1  —  t  “u,n  — 1  — 

9n-2  gn-2  9n-2  ffn-2 

=  jn  —  2Cv,n—2  ^  _  Un—2  jn  —  2Rv,n—2 

Cu,n  — 1  —  ?  -rtu^n  — 1  —  • 

9n-2  9n-2  9n~2 

The  unknown  Un-2  is  eliminated  from  Eu,n-\  by  writing 

Tn  — 1  —  *n  — l^u,n  — 2i  jn~l  *  Jn  —  1  "i"  in— l^u,n— 2> 

On  — 1  ~  in  — l^u,n  — 2)  Un  —  1  *  ^n  — 1  ^n  —  \Ru,n—2 

under  which  Eu,n-\  now  becomes 

9n  —  \'^n  "b  ^n  — l^n  — 1  d"  7n  — l^n— 2  d"  Jn  — l^n  — 1  d"  Otji  —  iVn  —  Un-l- 


Unknowns  Un-i  and  Un-2  nre  eliminated  from  this  equation  by  using  the  transforma¬ 
tions  _  _ 

jn  —  l  *  jn  —  \  d"  'Vn  — lOti.n— 2  d"  ^n  — inu,n  — 1) 

On  — 1  *  On  — 1  d"  T'n- 1  ^w,n— 2  d"  ^n- l^u,n  — Ij 

f^n  —  l  ~  d-'yn  — 1  Cv,n— 2  d"  iin  — l^u,n  — It 

Un  —  l  *  U n  —  l  d"  "fn  —  lRv,n—2  d"  hn  —  iRu,n  —  l' 

It  only  remains  to  eliminate  u„_2  from  E„^n,  using 


dfi  >  dti  -|-  enUn,n— 2>  Cn  ^  Cn  d"  ^n^v,n  —  2i 

bn  *  bn  d*  Cn^D^n- 2>  Un  ^  Un  ^nRv,n  —  2 

and  this  step  is  repeated  for  n  =  5, 6, . . . ,  iV.  Employing  the  upper  boundary  conditions 
gives  VN+i,  vn  and  un-h,  and  then  the  modified  Eu,n-i  implies  that 

Un-i 

VN-l  =  - - ■ 

jN-1 

Using  (i46a)  with  n  =  N  gives  vyv-2  =  (un-i> I'/v-l-i)  and  (^66)  yields  u;v-i  = 
<t>2  (vN-ii  Viv+i)-  In  this  way  it  is  straightforward  to  back  substitute  and  find  the 
solution. 


33 


APPENDIX  B 


In  this  appendix  we  summarise  the  algorithm  used  to  solve  the  inhomogeneous 
equations  given  in  (3.5).  The  method  mimics  the  algorithm  for  the  homogeneous  case 
but  with  slight  adaptation  to  allow  for  different  boundary  conditions.  Additionally, 
the  resulting  system  can  be  solved  directly  and  does  not  require  any  iterations. 

The  boundary  conditions  at  y  =  0  (i.e.  u  =  v  =  v'  =  0)  are  used  to  initialise  the 
technique;  first  it  is  noted  that 


92  ,  U2 

U2  — 

h.2  /12 


It  is  then  possible  to  eliminate  this  quantity  from  the  equation 
formations, 


U3 


using  the  trans- 


The  next  equations  may  now  be  solved  to  obtain  the  relationships 


V3  =  av,3V4  +  6t,,3U5  +  Ct,,3U4  +  Rv,3,  (Blo) 

U3  =  Su,3V4  +  &tt,3V5  +  Sti,3t^4  +  Ru,3-  (-516) 

These  functions  are  eliminated  from  the  equations  £^,4  and  £^v,4  by  implementing  the 
transformations 


C4  — »  C4  -f-  d4av,3i  64  — ♦  64  +  d4  6v,3>  fi  fi  +  d4C„,3,  V4  — >  V4  —  d4Ry  3^ 

H  j\  +  HO‘U,3i  04  =  H^u,3i  /J4  — +  /14  +  140^,3,  U4  —*  U4  —  i4Ru,3- 

The  current  equations  may  now  be  solved  to  obtain  expressions  for  V4  and  U4  similar 
to  those  in  (Bl).  This  procedure  can  be  generalised  and  represented  as  follows.  Once 
expressions  similar  to  (Bla,b)  have  been  derived  for  Un-i,  Vn-i  and  Vn-2  then  write 

jn  ^  jn  ■f'  inO‘v,n  —  lt  Oln  =  inbv,n  —  li 

hfi  *  hfi  in(^v,n  —  l  7  Un  ^  bJji  InRv ,n  —  l 

and 

dfi  >  dfi  -}-  Cft  ^  Cji  -|-  Bnbv,n  —  27 

In  ~  ^n^v,n  —  27  In  *  In  ^nRv,n  —  2 
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where  7„  is  the  coefficient  of  u„_i  in  To  eliminate  u„_i  and  u„_i  from  Ev,n  the 

transformations  are  used 

Cjv  *  Cn  d"  dnOv.n  — 1  d"  7n®u,n  — Ij  ^  d"  dn6u,n  — 1  d"  7n^ti,n  — 1? 

/^n  ^  /"n  d“  dnCv,n  —  l  d"  ')n^u,n  —  \t  '^n  *  —  i  7n-Ru,n  — 1) 

are  used.  The  forms  of  the  double  over-barred  terms  may  be  generalised  according  to 


^i;,n  — 


U  1  fn  Qn 

^  r  _ 

p  ^  —  p  ) 

n  A  n 


Rv,n  — 


1/  fnUn 

fen 

r  ’ 


=  _  tin  Jn^v,n 

~  ~T  ’ 

fin  hn 


^u.n  — 


fnSn 

Tnhn 

jn^v,n 


hn  ’ 


=  9n  jn(^v,n  "p  _  jn^n 

Cu,n  —  ~T  r  >  iiu,n  —  T  7  ) 

fin  fin  iln 


with  r„  =  C„  -  fnjnlhn. 

This  procedure  is  repeated  until  the  upper  limit  of  the  solution  space  is  reached 
(entirely  analogous  to  the  method  given  in  Appendix  A).  It  is  then  elementary  to  back 
substitute  and  deduce  the  solution. 

In  the  solution  of  the  small  wavenumber  equations  (6.2)  it  was  known  from  the 
results  of  Hall  (1985)  that  the  most  dangerous  mode  for  this  system  corresponds  to 
eigenfunctions  which  are  even  in  the  co-ordinate  ip.  In  order  to  solve  (6.2)  it  was  there¬ 
fore  convenient  to  only  consider  the  domain  [0,  oo)  and  alter  the  boundary  conditions 
at  ^  =  0  in  order  to  ensure  that  an  even-valued  eigenfunction  set  was  indeed  derived. 

In  the  process  of  doing  this  it  was  important  that  the  third  derivative  of  a  function 
was  calculated  accurately  when  the  iteration  process  on  the  desired  eigenvalues  was 
applied.  The  discretisations  discussed  in  Appendix  A  were  found  to  be  insufficiently 
precise  for  this  purpose  and  therefore  higher  order  discretisations  were  required.  For 
example,  the  second  derivative  of  a  function  at  the  origin  was  now  written  as 
(P <p  2,(f)\{h\h2h3  -|-  hih^h^  -|-  /i2/i3/i4  -|-  -|-  h\hnh^  -|-  /i2/j4/i5  -|-  h^hih^) 


dy'^  h\h2hzhih^ 

2^2(^2^3^4  +  ^2^3^5  +  /l2^4^5  +  h^h^hp) 
h\{h2  —  h\){h3  —  /ii)(/i4  —  ^i)(/i5  ~  hi) 
^  2(^3 (/i  1/13/14  -f-  hih^h^  -t-  hih^h^  -f-  h^hih^) 
[hi  —  /l2)/*2(/*3  ~  h2){hi  —  /l2)(/i2  ~  /is) 
2(^4(/ll /12 /i4  "k  /il/l2/i5  "t"  hih^h^  -f-  /l2/j4/i5) 
{hi  —  /i3)(/i2  ~  /i3)/i3(/i4  ~  /i3)(/i3  —  /^s) 
^  2(^5(/i1 /l2/l3  +  /ll/i2/i5  +  hih^h^  -t-  /l2/i3/i5) 
{hi  —  hP){h2  —  /i4)(/i4  ~  /i3)/i4(/i5  ~  h\) 
2(^6(/il /i2/i3  "b  /ii/i2/i4  +  hih^hi  -j-  /l2/l3/i4) 
{hi  —  /l5)(/i2  ~  /i5)(/i3  ~  /i5)/i5(/i5  ~  h^) 
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where  hj  =  yj^i—yi .  Similar  expressions  were  used  for  the  third  and  fourth  derivatives 
at  the  origin  and  neighbouring  points  but  they  are  not  included  here  for  reasons  of 
brevity.  However  details  may  be  obtained  upon  application  to  either  author. 
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(o)  Oo=3,  c, 


(b)  Qo  =  2,  c 


Figure  3.  The  coefficient  cj  of  amplitude  equation  (3.9)  as  a  function  of  the  vortex 
wavenumber  k  for  ffo  =  (a)  3,  (6)  2,  (c)  1,  (d)  0,  (e)  —1,  (/)  —1.5  amd  (g)  —2. 


(a)  no  =  3,  C2 


(b)  Qo  =  2.  C2 


Figure  4.  The  coefficient  C2  of  amplitude  equation  (3.9)  as  a  function  of  the  vortex 
wavenumber  k  for  fio  =  (<*)  3,  (6)  2,  (c)  1,  (d)  0,  (e)  —1,  (/)  —1.5  and  {g)  —2. 
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Figure  4.  The  coefficient  cj  of  amplitude  equation  (3.9)  as  a  function  of  the  vortex 
wavenumber  A:  for  Qq  =  (a)  3,  (b)  2,  (c)  1,  (cl)  0,  (e)  —1,  (/)  —1.5  and  (ff)  —2. 


Figure  5.  The  equilibrium  amplitude  ’I^<1  (defined  by  (5.1))  as  a  function  of  the 
vortex  wavenumber  k  for  Qo  =  (a)  3,  (6)  2,  (c)  1,  (d)  0,  (e)  -1,  (/)  -1.5  and  {g)  -2. 
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Figure  6.  Schematic  diagram  of  the  asymptotic  structure  of  the  solution  of  homo¬ 
geneous  system  (3.3)  for  the  case  of  small  vortex  wavenumber  A’  <C  1.  The  solution 
configuration  divides  into  four  distinct  zones:  Zone  /  is  a  thin  region  of  depth  0(A^ ) 
at  a  distance  ~  1.47fio  ^’2  from  the  wall.  This  zone  is  embedded  within  //,  a  region  of 
thickness  ©(As).  The  structure  has  an  additional  viscous  layer  at  the  wall,  zone  III, 
which  enables  the  boundary  conditions  there  to  be  met.  Finally,  a  region  IV  of  depth 
0(1)  facilitates  the  exponential  decay  of  the  disturbance  solutions  far  from  the  wall. 
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